Simulating high-temperature superconductivity model Hamiltonians with atoms in 

optical lattices 
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We investigate the feasibility of simulating different model Hamiltonians used in high-temperature 
superconductivity. We briefly discuss the most common models and then focus on the simulation of 
the so-called t-J-U Hamiltonian using ultra-cold atoms in optical lattices. For this purpose, previous 
simulation schemes to realize the spin interaction term J are extended. We especially overcome 
the condition of a filling factor of exactly one, which otherwise would restrict the phase of the 
simulated system to a Mott-insulator. Using ultra-cold atoms in optical lattices allows simulation 
of the discussed models for a very wide range of parameters. The time needed to simulate the 
Hamiltonian is estimated and the accuracy of the simulation process is numerically investigated for 
small systems. 

PACS numbers: 03.75.Ss, 03.67.-a 



I. INTRODUCTION 

Since its discovery, high-temperature superconductiv- 
ity has attracted much attention on the theoretical 
as well as on the experimental side. Nevertheless, the 
mechanisms which lead to this effect are not completely 
understood Q{. On the theoretical side one encounters 
the problem that the physics of a very complex many- 
body system has to be described and calculated. Because 
quantum effects play a crucial role, the calculation of 
these systems is a very hard task on a classical computer. 
Furthermore it is not certain that the investigated the- 
oretical models accurately describe the macroscopic fea- 
tures of superconducting materials used in experiments. 
Experimental efforts to test these models have suffered 
from difficulties in changing the system parameters for 
a given superconductor. In general different parameter 
sets require different types of superconducting material. 
Due to all these obstacles it is worthwhile searching for a 
versatile system which can be used to accurately simulate 
high-temperature superconductor models. 

As we will show in this paper ultra-cold atoms in op- 
tical lattices are a very good candidate for simulating 
high-temperature superconductivity. During the past 
few years the field of cold atom physics has made tremen- 
dous progress and entered the regime of accurately con- 
trolled strongly correlated systems. Seminal experiments 
with bosonic atoms have, e.g., realized the Mott-insulator 
to superfluid transition and cold controlled collisions be- 
tween atoms in optical lattices 0, 0, IE 01 • They found 
excellent agreement with corresponding theoretical mod- 
els 0, • For clouds of fermionic atoms superfluidity has 
been demonstrated by creating vortices [j| and Fermi sur- 
faces of atoms loaded into three-dimensional optical lat- 
tices have been measured |l0|. Optical lattice systems 
are also very flexible; a large range of Hamiltonians can 
be realized and system parameters are easily varied over 
a wide range by changing external laser parameters . 
For instance, it has also be en p roposed to implement ef- 
fective magnetic fields El El El El El E3 and even 



the implementation of non-Abelian gauge fields seems to 
be feasible E1E3 

We will briefly discuss some of the model Hamiltonians 
which are used to describe high-temperature supercon- 
ductive cuprates. We will then focus on the t-J-U Hamil- 
tonian 20], which includes a spin-spin coupling term J 
additionally to hopping t and on-site interaction U . The 
simulation of this Hamiltonian can be used for obtaining 
a deeper understanding of the corresponding model, and 
also enlarges the parameter space by an additional, freely 
tunable interaction term J. In analogy to the cuprates we 
expect this model to show a very rich phase diagram for 
the atoms in the optical lattice including the appearance 
of (anti)ferromagnetic phases. 

The simulation of the spin-spin interaction part in op- 
tical lattices has already been discussed in Refs. pll I22 . 
l23j | . Atoms with two internal states simulate a spin chain 
and the spin-spin interaction is implemented by state de- 
pendent shifting of the atoms causing controlled atomic 
collisions. These collisions induce state dependent phase 
shifts of the atomic wave function mimickin g^ th e spin- 
spin interaction. However, all of the proposals |2lLl22ll2^ | 
require a filling of very close to or exactly one atom per 
lattice site. For superconductors, this corresponds to a 
Mott-insulating state. In order to extend this method to 
the simulation of the t-J-U Hamiltonian we will relax the 
condition of having exactly one atom in each lattice site. 
Additional shifts of the atoms will be necessary for simu- 
lating the spin interaction term J and we will also show 
that by choosing appropriate parameters it is possible to 
implement effectively attractive interaction terms U < 0. 

The simulation of the time evolution of the whole t- 
J-U Hamiltonian is performed via a Trotter-Suzuki ex- 
pansion H3, E3, EE H3, H^l and we will consider the 
case of temperature T = in our calculations. To be 
able to experimentally observe superconductivity effects 
using our simulation method the temperature of the op- 
tical lattice atoms will have to be lower than IcbT < 0.02i 
|29| . where fee is the Boltzmann constant. Such low tem- 
peratures are experimentally difficult to achieve. How- 
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ever, theoretical proposals for fault tolerant loading of 
fermionic atoms |3f| into the lowest motional band of 
an optical lattice and phonon assisted side band cooling 
within this motional band [3l| exist. These methods will 
be realizable with current and near future experimental 
techniques and enable achievement and control of tem- 
peratures sufficiently small for our purpose. 

This work is organized as follows. In Sec.|H]we briefly 
discuss some of the model Hamiltonians used to de- 
scribe high-temperature superconductivity in cuprates. 
In Sec. IIIII we explain how one of these Hamiltonians, 
namely the t-J-U model, can be simulated using ultra- 
cold atoms in optical lattices. The simulation process is 
numerically tested in Sec. II VI and we conclude in Sec. El 



II. MODEL HAMILTONIANS FOR HIGH-T C 
SUPERCONDUCTORS 

In this section we briefly discuss some of the 
model Hamiltonians which are used to describe high- 
temperature superconductors [3^, 03- The most 
basic microscopic description is given by the Hubbard 
model 

H tu =~t ~A,<j~ c i-<y + U Y1 n^h jd . (1) 

The operator ej a (£j. a ) creates (annihilates) an electron 
with spin a =|, J, in lattice site j and obeys the stan- 
dard fermionic anticommutation relations, the number 
operator is denoted by hj l<7 — c^ a Cj^ a . This Hamilto- 
nian describes electrons in the lowest Bloch band tightly 
bound to the lattice sites formed by the atoms of the 
solid. The electrons can tunnel from one lattice site to 
the nearest-neighbor site (indicated by the brackets (i, j)) 
with hopping energy t and the repulsive on-site Coulomb 
interaction is given by U . Because of the Pauli exclu- 
sion principle for electrons and the single Bloch band 
assumption in Hamiltonian Eq. there can be a max- 
imum of one electron of each spin state in a single lattice 
site. Hence no onsite interaction term for electrons in the 
same spin state is included in the Hamiltonian. 

The Hamiltonian Eq. Q describes the behavior of the 
electron gas in two-dimensional (2D) layers. If the lattice 
sites are all occupied by a single electron the system be- 
haves like a Mott-insulating anti-ferromagnet. By hole- 
doping of the cuprate the behavior of the gas changes con- 
siderably and for a critical doping and temperatures be- 
low T c the solid gets superconductive (see, e.g., [29ll35|). 
In a real superconductor many of these two-dimensional 
layers are stacked on top of each other. The influence of 
these layers can be included by a small inter-layer tun- 
neling strength |3t| . 

Quantum cluster calculations give very strong evidence 
that the Hubbard Hamiltonian Eq. Q is already a suf- 
ficient minimal model to describe the main properties of 
the phase diagrams of cuprates [37j • Nevertheless, these 



calculations suffer from finite size effects. The simulation 
of this Hamiltonian using ultra-cold fermionic atoms in 
an optical lattice is straightforward and can overcome 
this problem for sufficiently large lattices. It requires the 
loading of the atoms into the lattice and choosing the 
lattice parameters such that two-dimensional layers are 
created. Details on which parameter range can be simu- 
lated are discussed in Sec. IIII Al 

However, the Hubbard Hamiltonian Eq. does not 
describe all properties encountered in a superconduc- 
tor [33. Other effective model Hamiltonians have been 
discussed trying to incorporate such experimentally ob- 
served effects. For example, it has been suggested [3!| to 
introduce a next-nearest neighbor hopping term 

H nn = -t' J2 Ca- (2) 

These terms can be realized in the simulation process by 
using very shallow optical lattices. In this case, however, 
the probability to excite atoms into higher Bloch bands 
increases and off-site interaction terms U os are important. 
Also it becomes more difficult to adjust all parameters in 
this Hamiltonian independently. Therefore, in this pub- 
lication, we do not discuss the realization of this model 
Hamiltonian any further. 

Another effective model is the so-called t-J-U Hamil- 
tonian, which has recently been used in connection with 
gossamer superconductivity [2(1 ■ It is given by 

H = H tu + Hj , (3) 

where 

Hj = jJ2s i -S j (4) 

is the spin interaction Hamiltonian of electrons in neigh- 
boring lattice sites. The spin operator Sj is defined by 
its components 

S W = l E 5 ^C, (5) 

?>c=u 

where a w are Pauli matrices, w = x,y, z. For very large 
interaction strengths U Hamiltonian Eq. J3J converges 
to the so-called t-J model, where double-occupancy of a 
single lattice site is excluded. This Hamiltonian has also 
very often been used to describe the behavior of high-T c 
cuprates. In the following sections we will investigate in 
detail how the t-J-U Hamiltonian can be simulated using 
atoms in optical lattices. Such simulations can help to 
understand and test features of the model, and also to 
investigate a larger parameter space which is likely to 
exhibit very interesting phase diagrams. 

III. SIMULATION OF THE t-J-U MODEL WITH 
ATOMS IN OPTICAL LATTICES 

In this section we show how Hamiltonian Eq. can 
be simulated by using ultra-cold atoms in optical lattices. 
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First we explain how to implement the [/-Hamiltonian 
and the spin-interaction independently from each other. 
Then these Hamiltonians are combined to realize the full 
t-J-U Hamiltonian. 



A. The hopping and interaction terms 

The Hamiltonian H t u, Eq. is simulated using a 
three-dimensional optical lattice where hopping alon g th e 
z-direction is suppressed by high potential barriers 
The lattice is filled with ultra-cold fermionic atoms occu- 
pying the lowest motional Bloch band only and moving 
in planes parallel to the sy-plane. We restrict our consid- 
erations to one such 2D layer of the optical lattice. The 
two spin states of the electrons are represented by two 
long-living internal states of the atoms which we denote 
by a =T,|- The Hamiltonian which describes the dy- 
namics of fermionic atoms in the lattice is given by Htu, , 
cf. Eq. QJ, with the hopping constant t and the interac- 
tion strength U s . By an appropriate choice of the lattice 
constants the system parameters t and U s can be tuned 
over a wide range and using a Feshbach resonance 
it is possible to change the s-wave scattering length be- 
tween the atoms which gives additional control over the 
interaction strength U s . 

However, we are not completely free in our choice of 
the lattice depths. In order to observe the superconduct- 
ing phase in high-T c cuprates the temperature has to be 
lower than the critical temperature T c , which is for high- 
temperature superconductors of the order of 0.02i |29| . 
For our system with atoms in an optical lattice we expect 
the same behavior. We consider a lattice depth of 5-Er 
with En = h 2 (2ir) 2 /2mA 2 the recoil energy, m the mass 
of the atoms and A the wave length of the laser produc- 
ing the lattice potential. In this case the hopping term 
t is of the order of 0.07-Er and thus the superconduct- 
ing phase will be observable for temperatures lower than 
approximately 0.0014£"ri/fcB- For the fermionic species 
Li trapped in a lattice with A = 670nm this corresponds 
to a temperature of 5nK. In recent experiments with 6 Li 
temperatures of about 30nK have been reported [j| . The 
required lower temperatures can be reached for example 
with phonon side band cooling [30l l3l| . 

As an aside we note that in the case of bosonic atoms 
we have to assume a very high interaction strength 
f/p = U[[ between two atoms in the same state to re- 
alize the above Hamiltonian H t u s and replace anticom- 
mutator relations with commutator relations. This in- 
teraction has to be much larger than the inter-species 
interaction strength [/-p = U s , the hopping constant t 
and the spin-spin interaction strength J. In this case all 
states with two or more identical atoms in the same lat- 
tice sites can be discarded. However, a direct mapping 
to a fermionic system using a Jordan- Wigner transforma- 
tion is only possible in one spatial dimension [41| . There 
are proposals for similar mappings in higher dimensions 
which require the addition of Majorana fermions to the 



system 0> 03 ■ Although appealing from a theoretical 
point of view it is not clear how to realize such systems 
with atoms in optical lattices. Nevertheless, in the case of 
bosons the above conditions on the interaction strength 
lead to a rich phase diagram similar to models discussed 
in Ref. H|. 

B. Implementation of the spin-spin interaction 

In general, the van der Waals interaction between two 
atoms in neighboring lattice sites is not sufficient to re- 
alize the required ratios of J, t, and U . Lowering the 
lattice barriers such that there is a significant nearest 
neighbor interaction also causes problems since in this 
case it is more likely that higher Bloch bands are occu- 
pied. Hence the spin interaction J has to be simulated 
indirectly. Methods for achieving this have already been 
detailed in Refs. |2lJ, |22j, |23J . These proposals require a 
filling factor of very close to or exactly one atom per lat- 
tice site which corresponds to half-filling of the electron 
system and thus a Mott-insulating phase |35|. Therefore, 
these proposals cannot be used to simulate superconduc- 
tivity. 

In order to circumvent this problem we extend the 
scheme proposed in Ref. [22]. Let us first consider the 
one-dimensional spin-spin interaction in z-direction, i.e., 

H zz — J z ^ Sj 

j (6) 

{i,3) 

This type of interaction between the atoms is realized by 
state selectively moving atoms. They are stored in very 
deep optical lattices such that any hopping is strongly 
suppressed, i.e., t = 0. The atoms in, e.g., state | f) 
are shifted and overlapped for a certain time with their 
neighboring atoms to the left and to the right in state 
| J,). The resulting atom- atom interaction leads to the 
desired phase shift as described in Ref. For illustra- 
tion let us assume two lattice sites described by the Fock 
states | nip nj; np n 2 ), where n J a is the number of atoms in 
state \o) in the jth lattice site. After the above shifting 
procedure this state evolves to 

|npnj;npn|) — > exp(ix(n}n 2 + n}n 2 ))|nj, n\; np n 2 ) , 

(7) 

where x is the phase acquired during the collision pro- 
cess. This is exactly the action of the last two summands 
—n^h J , — n\nl of Eq. JBJ. 

However, this process is not yet sufficient to imple- 
ment the spin interaction in z-direction. In order to 
achieve our goal the spin state of the atoms in every 
second lattice site has to be flipped, i.e., the operation 
Vfl = of ® I2 ® erf ® I4 <g) ... has to be applied to the 
atoms. This operation requires addressing each second 
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lattice site which can be done by using an additional 
standing wave laser field. The resulting supcrlattice must 
have twice the wavelength of the original lattice and its 
intensity minima need to coincide with every second lat- 
tice site of the trapping potential |44| . Thus the energy 
levels of every second atom are AC-Stark shifted out of 
resonance such that a microwave or laser-field driving the 
transition | f) | J.) realizes Vs.. 

By repeating this shifting process after the spin flip op- 
eration Vs and flipping the atoms back the state evolves 
according to 

K' n l; n T' n l) — * cx P(- i x(«l«i+«T n T))l"'T' n l' n T' n l) > 

(8) 

where the collision time has to be chosen such that the 
acquired phase is equal to — x- The whole process induces 
dynamics according to Hamiltonian Eq. ©, where x = 
J z r/4ft and r is the time for which H zz is applied (details 
on calculating the phases can be found in Appendix |A"|) . 

The creation of such spin-dependent phases has al- 
ready been demonstrated experimentally. In Ref. Q 
the authors used Rb atoms in an optical lattice with 
Vo = 34_Er. Excitations of atoms into higher Bloch bands 
were avoided by using an appropriate time of r s = AOfis 
to shift the atoms spin-dependently into the neighboring 
lattice sites 6J. After holding the atoms in the shifted 
position for Th = 450/iS and shifting them back again 
collisional phases of x ~ 27r were achieved. 

With being able to implement Hamiltonian Eq. JJjJ it 
is also possible to realize the Hamiltonians 

H xx = J* ^ > w 

{i,3) 

Hyy = JyJ2^Sl ( 10 ) 

We first observe that with suitable laser pulses it is possi- 
ble to implement the following rotations on single atoms 

V y = exp (z V) , (11) 
V x = exp {-i\(T x ) ■ (12) 

By simultaneously applying one of these rotations on all 
atoms, implementing Hamiltonian Eq. JfJJl and applying 
the Hermitian conj ugat e of the same rotation we get the 
time evolutions |2lL 1221 l4q 

Vj exp V y = exp (-iSxx^j , (13) 

V% exp V x = exp • (14) 

The interaction strengths J x , J y , and J z can be tuned 
independently from each other by the choice of the col- 
lisional phase x during the simulation of the respective 
Hamiltonian. 



C. Combining the Hamiltonians 

We simulate Hamiltonian Eq . (f51) usin g th e well known 
Trotter-Suzuki expansion |H |3 H Hi El E| de- 
scribed in detail in Appendix [5] In this approach the 
different parts of the Hamiltonian are simulated for a 
small time r separately. For instance, in first order, the 
decomposition of the t-J-U Hamiltonian reads 

° XP ( Z f I ) = 6XP {^hj ^ 6XP i^ 1 ) V ° 

x VJ exp (=^j V y exp (=^§^j + ^(r 2 ) . 

(15) 

The last part exp C—iHtu, T /hj of this time evolution 
is simulated by choosing a suitable lattice depth such 
that the required values of t and U s are realized. After 
waiting a time r |4fij the lattice has to be ramped up 
avoiding any excitation into higher Bloch bands. Then 
the shifting and flipping processes Vs, as described in the 
previous section, have to be applied in order to simulate 
the time evolution of H zz , H yy , and H xx . The latter two 
require the implementation of the operations V x and V y , 
cf. Eqs. (|llfl and 1)12(1. By repeating all of these sim- 
ulation steps m times a time mr is simulated with an 
error oc r 2 m. Better accuracies can be achieved by using 
higher order approximations, see Appendix iBl 

During the simulation of the spin-spin interaction term 
two atoms may occupy the same lattice site j. This gives 
rise to an additional term not considered so far. We 
illustrate this by considering shifts in ^-direction during 
the simulation of the spin interaction. The atoms will 
acquire an additional phase (compare Appendix 

f T0 ( K(r)-^(r)] 2 \ 

4> = Kj^dre X p(- [ TV J 2x , lK j , (16) 

with K given in Eq. I|A3(1 . This phase corresponds to an 
(additional) effective interaction term 

Hf = U eS J2n^h jtl , (17) 
j 

where the interaction constant is given by C/ e ff = tkp/r. 
Note that this effective interaction does not diverge for 
r — * as in this case also <j> goes to zero. This interaction 
provides an additional opportunity for tuning the simu- 
lated interaction constant U . The total interaction for 
small U c ft <C 2ttH/t is given by U — U s + 3dU c s with d 
the number of spatial dimensions in which the spin-spin 
interaction is simulated. Here the factor of 3 in front 
of the dimension number d arises because the phase <f> 
occurs in the simulations of H xx , H yy , and H zz which 
are all simulated for the same time r. By an appropri- 
ate choice of U s and 4> it is possible to simulate effective 
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FIG. 1: Lower bounds for the calculated shifting and holding 
times T s h (cf. Eq. 1181 ') for lattices with 87 Rb atoms (A = 
826nm, a s — 5.1nm), 40 K atoms (A = 826nm, a s = 5.5nm) 
and 6 Li atoms (A = 670nm, a s — 2.4nm). The values given 
in brackets are typical values used for the calculation of r s h ■ 



interactions which are attractive. For this purpose, the 
phase <f> has to be chosen in such a way that it can be 
written as (j) — 2tt + U^t/H with a small negative U' eS 
and 2tt ^> \U^ s \t/H. The total simulated interaction is 
then given by U = U s + 3dU^ s . For not too large U s this 
leads to a negative U |47|. 

The total time for which the system can be simulated is 
restricted by the life- and decoherence-time of the ultra- 
cold atoms in the optical lattice. This is typically on 
the order of one second 0. Hence we have to estimate 
the time which is needed for simulating one time step of 
Hamiltonian Eq. J3J). During the shifting process heating 
of the atoms has to be avoided. Since the lattice is re- 
quired to be very deep during the atoms are shifted hop- 
ping is strongly suppressed. This leads to a flat Bloch 
band and thus no excitation can take place within the 
band during the shift. Any heating thus means that 
atoms are excited into higher Bloch bands. In order to 
avoid such excitations we require for each shift of atoms 
in, e.g., state | f) into the neighboring lattice site a time 
larger than the inverse of the lattice site trapping fre- 
quency Lu t = \/Vo/Ek h{2ir) 2 / m\ 2 as already demon- 
strated experimentally in Refs. |<5, 6]. Furthermore, the 
atoms have to be held in the shifted position for a certain 
time in order to achieve the desired phase. This phase is 
either x or 2n — x- Thus the average time for one shifting 
and holding process is given by 
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The minimal shifting times for varying lattice depths are 
shown in Fig. ^ and are of the order of a few hundred 
microseconds for Rb or a few tens of microseconds for Li. 
To simulate the time evolution of H zz these shifts have 
to be repeated 2d times. 

The constraints for the ramping time between the ap- 
plication of the i-[/-Hamiltonian and the spin-spin inter- 
action are twofold: The ramping has to be adiabatic on a 



time scale given by 1 /ui t to avoid excitations into higher 
Bloch bands and it has to be rapid compared to Hn/2t. 
This ensures that the system does not evolve for too long 
with a hopping term t that is different from the desired 
one or even follows the change in t adiabatically in con- 
trast to the sudden change required by the Trotter-Suzuki 
expansion. The first constraint has been experimentally 
tested for Rb atoms |4j and was found to be fulfilled 
for ramping times longer than 50/xs in accordance with 
theoretical calculations |4q|. For typical lattice param- 
eters the hopping time is given by tvn/2t = 1ms and is 
thus more than one order of magnitude larger than the 
minimal ramping time required to avoid excitations into 
higher Bloch bands. Therefore, both constraints can be 
fulfilled at the same time. Furthermore, the ramping 
process conserves quasi-momentum, and excitations that 
do not change the quasi-momentum involve at least two 
hopping processes. Because the lattice is ramped on a 
time scale short compared to tm/2t excitations within 
the lowest Bloch band are, therefore, also strongly sup- 
pressed. Together with the time needed to simulate the 
t-J7-Hamiltonian it is thus feasible to simulate several 
hundreds of Trotter steps m within the lifetime of the 
atoms in the lattice. 

We finally remark that this procedure for simulating 
the t-J-U Hamiltonian is compatible with the additional 
simulation of a magnetic field. Various methods for cre- 
ating effective magnetic fields at the same time as the 
hopping and interaction terms in H t u were recently pro- 
posed miii ehii mill 

Thus it is possible to extend 
our setup for studying the t-J-U Hamiltonian in external 
magnetic fields. 



D. Measuring the properties of the atom gas 



The properties of the time-evolved state can be probed 
by measurements of first- and second-order correlation 
functions as proposed in Ref. 0. For this purpose 
a multiple matter-wave version of the Hanbury Brown 
and Twiss experiment [50l | has to be realized. The 
atoms are released from their trapping potential and 
imaged after a certain time of flight. Depending on 
the phase of the atomic gas these images reveal inter- 
ference patterns; second-order correlation functions give 
additional information on the phase and the underlying 
structure of the lattice. Such measurements have already 
been performed for bosonic atoms in the Mott-insulating 
phase 0, l5l[ and agree very well with the theoretical 
predictions (52}. With these techniques it will be possi- 
ble to check the phase diagrams of the simulated systems 
and compare them to calculated or measured phase di- 
agrams for different high-temperature superconductors 
(see, e.g., jaiH E3). 
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FIG. 2: Simulation of Hamiltonian Eq. (J^J using first- and 
second-order Trotter expansions, cf. Eqs. 1BI1 and 1B2II . Up- 
per bounds for the anti-fidelity T are shown. The number of 
lattice sites is M — 5, the parameters for the simulation are 
J = 0.3t, U s = 5t, U' cS = -2t, m = 1 leading to U = -t. 

IV. NUMERICAL SIMULATIONS 

As already mentioned in the introduction the simula- 
tion of Hamiltonian Eq. @ is a very hard computational 
task. Even for a few lattice sites an exact simulation 
of the system is no longer feasible and one has to use 
approximations such as mean-field theories or Quantum 
Monte Carlo calculations. Since we want to compare the 
results of the simulation of Hamiltonian Eq. J2J with ex- 
act results we have thus restricted ourselves to the one- 
dimensional case with a few lattice sites only. 

The anti-fidelity of the simulated state compared to 
a state which is time-evolved using the full Hamilto- 
nian Eq. is calculated as follows. Let Uq(t) = 
exp(— ii/r/fi.) be the full time evolution operator, U(t) 
the simulated time evolution, \ip)- m an input state and 
= U(T)\ip) in , |r/>o) = %(V)|^)in. Due to the Cauchy- 
Schwarz inequality and the properties of the matrix norm 
we get 

\\U{r) -Uo{rW >|«VI " (VblXM " (m 
= |2 - 2Ke«V#o»| • 1 1 

Since | (ip\ipo) | < 1 the rhs of Eq. I|19|l is always positive 
even without the modulus and we can rearrange 

\m Q )\>l-\\\U{r)-U Q {r)\\ 2 . (20) 

For ||W(t) — Uo(t)\\ 2 < 2 this yields an upper bound for 
the anti-fidelity T 

T =\-\m*)\ 2 <\\U{t)-U q {t)\\\ (21) 

which can be easily calculated from the time evolution 
operators. 

We have calculated this anti-fidelity J- for several cases. 
First, we investigate a parameter set which gives a nega- 
tive (attractive) interaction strength U and compare this 
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FIG. 3: Upper bounds for the anti-fidelity T for first- and 
second-order Trotter expansions, cf. Eq. <B4t . The simulated 
time is fixed to r = WOh/t, the number of lattice sites is 
M = 5. We have used J = 0.3t, U s = lOf, and = 0. 



to the exact time evolution according to Eq. 10. Upper 
bounds for the corresponding anti-fidelity T are shown 
in Fig. |21 It has been assumed that no errors such as 
unprecise creation of the phase shifts \ or excitation of 
higher Bloch bands occur during the shifting process. For 
small times the slopes of the curves for the first and sec- 
ond order approximations agree very well with the predic- 
tions from first and second order expansions, cf. Eqs. 1|B1|) 
and (|B2|) . Note that the slopes shown are increased by 
a factor of two because we take the square of the matrix 
norm in order to calculate the anti-fidelity T. Figure |3 
shows that for small times the simulation of the time evo- 
lution of Hamiltonian Eq. agrees very well with the 
full time evolution but the desired simulation times of 
order lOOh/t can only be reliably achieved by repeating 
this simulation process. 

The dependence of the upper bounds of the anti- 
fidelity T on the number m of Trotter steps for a fixed 
time r = lOOh/t is shown in Fig.[3]for values correspond- 
ing to typical parameters in high-temperature supercon- 
ductivity [54| . As expected the fidelity gets better for an 
increasing number of steps. In second order approxima- 
tion for m = 500 steps the anti-fidelity is already smaller 
than 10~ 3 . For the same accuracy using the first-order 
approximation m — 900 Trotter steps have to be used. 
This means in first-order approximation the lattice has 
to be ramped up and down 900 times and H zz has to be 
simulated 2700 times with 1800 applications of V x . y . In 
second-order each single step is more complicated, but by 
choosing the simulation process in a judicious way only 
500 rampings of the lattice are required. By optimizing 
the process (cf. Appendix only 2500 simulations of 
H zz are required with 1500 applications of V x>y . 

To get an impression on how the anti-fidelity T scales 
with the number of lattice sites M we calculated it for dif- 
ferent values of lattice sites M in a ID lattice. The results 
are shown in Fig.0|for a simulated time of r = O.Olh/t. 
The results indicate that the anti-fidelity T does not in- 
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FIG. 4: Upper bounds for the anti-fidelity T vs number of lat- 
tice sites M for first- (pluses) and second-order (stars) Trotter 
expansion. The simulated time was r = Q.Olh/t and we have 
chosen J = 0.3t, U s = Wt and = 0. 



are longer than the typical time scales associated with 
the dynamics of the t-J-U Hamiltonian. Hence the prop- 
erties of the Hamiltonian can be made visible using our 
simulation method. We also showed that by an appro- 
priate choice of the lattice parameters and shifting times 
it is possible to create attractive on-site interactions U. 
Also, our method is compatible with the simulation of 
magnetic fields in the lattice. 

Furthermore, we simulated the simulation process nu- 
merically for small systems and compared the results 
with exact calculations of the time evolution. For these 
simulations the results show very good agreement be- 
tween the simulated and the real time evolution. Anti- 
fidelities T < 10 -3 can be achieved. Because of these 
results and of the experimental progress we are confi- 
dent that the proposed method will help to deepen our 
knowledge and understanding of the t-J-U Hamiltonian 
model. 



crease quickly with the number of lattice sites. Therefore, 
it should be possible to experimentally simulate systems 
of realistic size. 

Since the time evolution of the full Hamiltonian Eq. © 
is only simulated approximately in our scheme, the state 
vector of the atoms after the time evolution will also only 
approximate the real state. However, we expect errors 
due to this effect to be small, since the necessary mea- 
surements to distinguish the phases involve only first- 
and second-order correlation functions. The influence 
of small imperfections of the simulated Hamiltonian on 
these functions will in general be smaller than the influ- 
ence on the state vector estimated by the upper bounds of 
the anti-fidelity T . Thus the deviations of the simulated 
phase diagrams is expected to be smaller than implied 
by the anti-fidelities. 



V. CONCLUSION 

In the present paper we have shown that it is feasible 
to simulate different model Hamiltonians which are used 
to describe high-temperature superconductive cuprates. 
The simulation of the most minimal model, the Hubbard 
model Eq. QJ, can already be expected to deepen our 
understanding of the phase diagrams encountered in the 
cuprates. Since this model does not describe all exper- 
imentally measured effects and to enlarge the accessible 
parameter space which can be investigated using atoms 
in optical lattices we have especially discussed how to 
simulate the t-J-U Hamiltonian Eq. @, which requires 
the inclusion of an additional spin-spin interaction term 
J. 

In order to make the simulation possible we have ex- 
tended earlier proposals to create an effective spin inter- 
action between atoms to the case where the lattice is not 
fully occupied. We found that near future technology will 
allow to simulate Hamiltonian Eq. © for times which 
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APPENDIX A: CALCULATION OF THE 
COLLISIONAL PHASES 

The exact phase shift due to the shift and interaction 
process can be calculated analytically jfj . If we assume 
deep lattice potentials of V x>y>z , we can approximate the 
wave functions of the atoms located in one lattice site by 
Gaussians. For illustration we assume an atom of state 
| 1) in lattice site j and one of state | J.) in lattice site 
j + 1. The shift shall occur in x-direction. Let x j ^(t), 

xj +1 (r) denote the time-dependent x-coordinates of the 
atoms. Because of the collision process the state of the 
two atoms will evolve to 



|l,0;0,l)^e-^|l,0;0,l) 
with the collisional phase given by 



2x1 



X — K dr exp ■ 

Here we have defined 



K 



4ira s h (V^Pr / V 



A 3 \ Ei 



(Al) 



(A2) 



(A3) 



where a s is the s-wave scattering length between the two 
atoms in states | f) and | J.), Vq is the depth of the lattice, 
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x = {/E R /V \/2ir and E R = ft ! (2tt) 2 / 'm\ 2 is the recoil 
energy, m the mass of the atoms and A the wavelength 
of the lasers creating the optical lattice. The time To is 
chosen such that the entire shift process is contained in 
the time interval from — tq to tq. For typical experimental 
parameters @ of Vq = 34£ , r and 87 Rb atoms with a s — 
5.1nm the prefactor in front of the integral has a value 
of roughly K — 19rad/ms such that after a collisional 
time of a few hundred microseconds phase shifts of the 
order of 2ir can be achieved. For Li atoms the prefactor 
has a value of the order of K = 230rad/ms and hence 
only roughly a tenth of the time is necessary to achieve 
similar phase shifts. The exact phase shifts depend on 
the details of the shifting process itself, i.e., on the exact 
form of the functions x 3 a (r) . 



APPENDIX B: THE TROTTER-SUZUKI 
EXPANSION 



where x — —ir/h and r is the simulation time. The 
operators Qi(x) and Q2{x) are approximations to the 
time evolution operator of i? ar b with 



exp(H alh x) = Q^x) + 0(\x\ 2 ) = Q 2 (x) + OQx\ 3 ) . (B3) 



Higher order approximations to the time evolution of 
-ff ar b exist, but they either require the simulation of neg- 
ative time which is difficult for the hopping Hamiltonian 
H t jj, or they include terms which make the simulation 
unstable |25j ]. Therefore, we simulate longer times using 
a generalized version of the Trotter formula [2^] with m 
simulation steps, reading 



For an arbitrary Hamiltonian i/ ar b 



=1 Hj the 



Trotter-Suzuki expansion 0, 0, HE HE l27U2S| is found 
by defining 



exp(i? arb a;) = Q 3 ) 



O 



X\3 + 1 

m'j 



(B4) 



Q*{x) 



Qi( x ) = n e ^v( xH j) > 

Hf-ix/2 HfX -Hf-ix/2 p Hix/2 



(Bl) When using the second-order approximation Eq. I|B2|) 
in this formula the experimental procedure can be opti- 
mized by combining subsequent simulations of H\ since 
> ( B2 ) exp(H 1 x/2)exp(H 1 x/2)=exp(H 1 x), 
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